Quasielastic magnetic scattering of neutrons at the systems with heavy fermions 



00 
OS 

o> 

< 

00 
(N 



i 

CO 

c3 



I 

C 

o 
o 



> 

m 
oo 

o 

00 
On 

c3 



S3 
O 
O 



X 
S3 



Andrei S. Mishchcnko 
RRC "Kurchatov Institute" 123182 Moscow, Russia 
(February 1, 2008) 

The theory of the quasielastic magnetic scattering of neutrons at the spin liquid with RVB (res- 
onance valence bonds) correlations is presented. Calculations demonstrate that the dependence of 
the scattering cross section on the energy transfer reproduces experimental shape of quasielastic 
peak in heavy fermion systems. It is shown that Fermi statistics of the spin liquid elementary exci- 
tations leads to the oscillations of the quasielastic scattering total cross section as a function of the 
momentum transfer. 
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1. Since it was suggested to treat the state with 
heavy fermions (HF) as the spin liquid (SL) with RVB 
type of spin-spin correlations the validity of RVB concep- 
tion is still the subject of considerable discussions. Al- 
though conception of RVB type SL appeared to be fruit- 
ful for quantitative description of the thermodynamic Q] 
and low energy spectral || properties of HF systems 
there is no unambiguous proof of existence of RVB cor- 
relations. The basic stumbling-block on the way to the 
SL identification is the lack of decisive experiment which 
can give possibility to accept or reject RVB conception. 

According to the RVB scenario the system of localized 
spins transforms at temperatures close to Kondo tem- 
perature Tk into the half filled band of magnetic exci- 
tations with the bandwidth T* ~ T K The basic 
fingerprint of the RVB correlations, which distinguishes 
SL state from the system of localized spins, is the Fermi 
statistics of the elementary excitations. Therefore, the 
study of consequences of the different statistical distri- 
butions which characterize SL and localized spin regime 
is a proper way to reveal some properties which are pe- 
culiar only to the highly correlated RVB state. 

In the present paper I suggest the theory of the 
quasielastic magnetic scattering of neutrons at the RVB 
type SL and show that change of the Boltzmann statis- 
tics of localized spins to the Fermi statistics of the SL 
excitations results in the oscillatory behavior of the to- 
tal quasielastic cross section as a function of momentum 
transfer. 

2. The basic model describing the HF state is the 
Anderson lattice Hamiltonian of the f-ions which are hy- 
bridized with the conduction electrons. The Coqblin - 
Schrieffer canonical transformation |||J|] eliminates the 
hybridization term and the lowest crystal field doublet 
(7 = ± can be treated in terms of exchange interaction 
^mm' between the spin states of different sites m = 1, £ (1 
is the number of the elementary cell, £ is the basis vector) 
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Here /^ CT (/mo-) are creation (annihilation) Fermi oper- 
ators of the spin state a at the site m which are subject 
to the constraint condition /i, ./ m(T = 1. 

In the mean field approximation the RVB state at low 
temperatures is characterized by nonzero intersite aver- 
ages (/mff/m'cr} and effective Hamiltonian of the SL low 
temperature states takes the form 
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which formally resembles tight binding approximation for 
band spectrum. The constants A mm i are determined by 
the values of exchange integrals I^m' an d anomalous av- 
erages ||. The chemical potential /i in the mean field 
approximation can be deduced from the global constraint 
for A-ion lattice A^ 1 £ m J2 a (fLJ™*) = 1- Therefore, 
the RVB SL in the mean field approximation can be con- 
sidered as a half filled band of Fermi quasiparticles with 
the width ~ T* . The eigenstates | A) of the Hamiltonian 
(Q) can be expressed in terms of eigenvectors in the 
form of superposition jq] 
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(Here ip(r — 1 — £) is the spatial wave function of localized 
f-electron; | a) is the spin component of the wave func- 
tion and k A is the wave vector of the state | A).) The 
occupation numbers of the states | A) with the energies 
£\ at temperature T are described by Fermi distribution 
n A - {l + exp[(e A -^/T]}" 1 . 

3. The standard expression for the scattering function 
JTo| | per one magnetic ion with the momentum transfer n 
and energy transfer Hlo can be expressed in terms of the 
matrix elements of operator Q = exp{iK(l + £)}J 
and Cartesian components of unit vectors k a = n a j \ k \ 
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(b is the magnetic scattering length |10| ; J is the operator 
of total magnetic moment). In case of the scattering 
at the localized crystal field states the statistic function 
is connected with the Boltzmann distribution and 
depends only on the occupation numbers of the initial 
states A. On the other hand in case of Fermi statistics 
of SL excitations the occupation of the final state is also 
important and statistic function takes the form 



^a,a' = n\(l - n\>). 



(5) 



Inserting the eigenstates (^) into expression (Q) and ne- 
glecting the intesite matrix elements of operator J be- 
tween highly localized /-electron states one gets the fol- 
lowing expression for the scattering function 
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Here B(k) — b 2 (gjF(n) /2) 2 where gj is the Lande factor 
and F(k) is the formfactor of the magnetic ion. 

4. First of all it have to be noted that dependence 
of the scattering function (^J) on the energy transfer fitu 
and temperature qualitatively describes the experimen- 
tally situation with the quasielastic peaks in HF systems 
p| . The quasielastic peakshape in HF systems is treated 
usually by use of the phcnomenological expression jfl| 
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Characteristic feature of HF systems is nonzero width 
F(T) at zero temperature. 

For quantitative estimation let us consider the simple 
case of symmetric spinon band with the constant density 
of states V(e) = 1 for e G [-T*, T*] and V{e) = out of 
the defined range. The chemical potential of this band 
is defined by condition of half filling and equal to zero. 
After integrating of the scattering function (^|) over the 
angles (which corresponds to the scattering at a poly- 
crystal) one gets the following expression 
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R(k) J deV(e + hw)n(e)(l-n(e + huj)). (10) 
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FIG. 1. The energy dependence of the scattering function 
of the spin liquid with the bandwidth T* — 2 for T = 0.1 
(thick solid line), T = 0.25 (thin solid line), T = 0.5 (dashed 
line) and T = 2 (dotted line). 



The momentum dependence is determined by the factor 
R(k) and will be considered in the next section. Now we 
concentrate on the dependence of the scattering function 
on the energy transfer which for the defined above density 
of states T>(e) has the following form 
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As may be seen from FigJIIthe energy and temperature 
dependence of expression ( |ll|) is in qualitative agreement 
with experimental measurements of quasielastic response 
in HF systems. Moreover, the energy dependence (|ll| ) is 
similar for Hlo < T* to the phenomenological function 
(^). It is instructive that generally accepted for the ex- 
perimental data treatment fit of the formula ^ to the 
scattering function ( pr| ) gives T(T = 0) ~ T*. It have 
to be noted that breaks of the energy dependence at 
hu> = and restriction of S s i (k, hu) by the energy range 
[— 2T*, 2T*] can be avoided by use of less singular density 
of states and involving the damping of SL excitations, re- 
spectively. 

5. The energy dependence of the quasielastic scatter- 
ing at the SL can not be used as a decisive argument in 
favor of RVB model applicability for HF systems because 
there are lot of models which give the same result (e.g. 
p^| , ^3[ ). However, much more bright consequence of the 
RVB model, which can serve as decisive test, is the de- 
pendence of the scattering function on the momentum 
transfer. 

To demonstrate this consequence of the Fermi statistic 
let us consider two-sublattice system in which intersub- 
lattice interaction dominates, i.e. -Ajtj/t' = 5\yT* /2 for 

£ 7^ £ and negligible for £ = (Numeric calculations 
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show that these approximations, although simplify the 
analytic estimations, do not influence the character of 
presented effects.) In the framework of this model the 
four-fold degeneracy of the states in each elementary cell 
is reduced and the set of eigenstates consists of two two- 
fold degenerate eigenstates with the energies ±T*/2. The 
chemical potential is fi — 0. Let us assume for dehnitness 
sake that the basis vector £ is parallel to z-axis and the 
distance between the magnetic ions is d. 

To calculate the momentum dependence it is conve- 
nient to introduce the function 



T* 



d(huj)S s i(K, Huj), 
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which corresponds to the integral value of the cross sec- 
tion with respect to the energy transfer. This function 
represents the total quasielastic scattering on the lowest 
crystal field level. 

In the final expression for the SL scattering function 
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where 
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the Fermi statistics results in the additional factor 

D(k, T) = ~ 1 1 - tanh 2 cos( K d) j , (15) 

which is equal to unity in case of Boltzmann distribution. 
The factor D(n,T) leads at low temperatures to the os- 
cillation of the scattering cross section with the period 
k p = 27r/(dcos[Kd]). Important property for experimen- 
tal check is the dependence of the oscillations period on 
the momentum transfer direction. 

The scattering function for polycrystalline samples 
SJ{?(k) can be calculated by integration of expression 
( |l3| ) over the angles. To demonstrate the consequences 
of Fermi statistic I separated out into the factor S tot (K) 
the standard scattering function for Boltzmann statistics 
which only source of the momentum transfer dependence 
is the square of formfactor F(k). It is seen from the 
following expression 



Qtot 



1 - tanh 2 
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that integration does not lead to the total suppression of 
oscillations. Here 
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FIG. 2. Dependence of the ratio («) = St? (k) / S tot (k) on 
the momentum transfer module n for SL bandwidth T* = 1 
and n = 0. Upper panel: in the perfect lattice T = (solid 
line), T = 0.7 (dashed line) and T = 1.0 (dotted line). Lower 
panel: at zero temperature for A = 1 (solid line), A = 0.6 
(dashed line) and A = 0.2 (dotted line). Insert: Scattering 
function taking into account the formfactor of Ce 3+ ion for 
the interionic distance d = 4A. 



V = 



(17) 



is the measure of the one site anisotropy of scattering. 
The character of the oscillations is determined by the 
factor sm(K,d) / (nd) and function 



$(/cd) 



sm(nd) — kc?cos(kc?) 



(18) 



For large enough momentum transfer nd 3> 1 the func- 
tion <J>(fttf) is much less than factor sin(«;rf) / (nd) . It is 
seen from the upper panel of Fig. ^ that the temperature 
does not significantly suppress the oscillations, at least 
in the mean field approximation. Even at temperature 
~ T* , which is enough to destroy SL state, the oscilla- 
tions amplitude is still noticeable. 

One should note that the period of driven by the Fermi 
statistic oscillations is determined by the distance be- 
tween interacting magnetic ions d and have nothing to 
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do with the peculiarities of SL Fermi surface. The char- 
acteristic distance between the magnetic ions in the HF 
compounds is d ~ 4A Therefore, in spite scattering 
function decrease ^"'(k) ~ F 2 (k) it is possible to ob- 
serve several oscillation periods (see insert in the upper 
panel of Fig. |^). 

6. However, the influence of the lattice imperfection 
on the suppression of the oscillations is much more sig- 
nificant. In the simplest case, when coherency is not 
destroyed by the absence of one of the magnetic ions, 
the lattice imperfection results in different environments 
of basis ions. The change of environment leads to the 
change of the wave function of the crystal field states jl4} . 
Different wave functions of basis ions lead, in turn, to in- 
equality of matrix elements {a | J a | a')\ ^ (<r \ J a | cr'} 2 
for sites 1 and 2. Although expression for scattering func- 
tion is derived for general case I present less cumbersome 
formula for the isotropic one-ion scattering case (x, y 
and z directions arc equivalent and Cartesian indices are 
omitted in the following). Nonequivalence of the basis 
ions can be characterized by the quantity 



U = 



£ (a | J | + (a | J | a') 2 



(19) 



Within this simple model the influence of the lattice im- 
perfection on the scattering function 



St? ( K ) = S tot ( K )\ jl- A tank 
is determined by the factor 

A 1 ~ U 



sin(Kcf) 
nd 



1 + U' 



(20) 



(21) 



which is equal to unity in the perfect lattice and tends 
to zero for the total lost of coherency. The lower panel 
of Fig. H demonstrates the influence of the lattice imper- 
fection on the suppression of the oscillations. 

7. In conclusion, it have to be noted that dependence 
of the spin liquid neutron scattering cross section on the 
energy transfer can not be distinguished from the scat- 
tering on the relaxing spin. However, the oscillatory de- 
pendence of the cross section on the momentum trans- 
fer is unique for spin liquid of RVB type. Besides, in 
case of soft enough crystal field splitting Acf ~ T* one 
have to include into the Hamiltonian (Q) the terms which 
are nondiagonal with respect to crystal field levels ||. 
This terms lead to the oscillations of the inelastic scat- 
tering intensity. The oscillations of the inelastic scatter- 
ing intensity were experimentally observed |l5[ in Kondo 
- semimetal CeNiSn and quantitatively explained in the 
framework of spin liquid conception However, since 
many characteristics of the periodic systems demonstrate 
the oscillatory behavior in neutron scattering response 



(e.g. the peak width oscillations in paramagnets JT6| , |l7[ ] ) , 
it is important to find the property which is connected 
with RVB correlations unambiguously. 

As it is shown in the present paper, the oscillations 
of the total quasielastic scattering cross section is direct 
consequence of the Fermi statistic of RVB elementary ex- 
citations. Therefore, an experimental observation of the 
predicted effect might serve as notable argument in favor 
of the validity of spin liquid conception for the descrip- 
tion of heavy fermion systems. 
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